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ABSTRACT 

We make a quantitative prediction for the detection rate of orphan gamma-ray burst (GRB) afterglows as a 
function of flux sensitivity in X-ray, optical, and radio wavebands, based on a recent model of collimated GRB af- 
terglows. We find that the orphan afterglow rate strongly depends on the opening angle of the jet (roughly oc f?^ 2 ), 
as expected from simple geometrical consideration, if the total jet energy is kept constant as suggested by recent 
studies. The relative beaming factor b K \, i.e., the ratio of all afterglow rate including orphans to those associated 
with observable prompt GRBs, could be as high as b K \ > 100 for searches deeper than R ~ 24, depending on 
afterglow parameters. To make the most plausible predictions, we average the model emission for ten sets of af- 
terglow parameters obtained through fits to ten well-observed, collimated GRB jets, weighted by the sky coverage 
of each jet. Our model expectations are consistent with the results (or constraints) obtained by all past searches. 
We estimate the number of orphan afterglows in the first 1500deg 2 field of the Sloan Digital Sky Survey (SDSS) 
to be about 0.2. The relative beaming factor b le \ is rapidly increasing with the search sensitivity: b re \ ~ 3 for the 
SDSS sensitivity to transient objects in the northern sky (R ~ 19), ~14 for the past high-z supernova searches 
(R ~ 23), and ^50 for the sensitivity of the Subaru Suprime-Cam (R ~ 26). Predictions are made for the current 
facilities and future projects in X-ray, optical, and radio bands. Among them, the southern-sky observation of the 
SDSS (sensitive to transients down to R ~ 23) could detect ^40 orphan afterglows during the five-year operation. 
Allen Telescope Array would find about 200 afterglows in a radio band at ^0.1-lmJy with b K \ ~ 15. 

Subject headings: gamma rays: bursts — ISM: jets and outflows — surveys 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are now confirmed to be located 
, at cosmological distances, and they are recognized as the most 
energetic explosion in the universe in terms of isotropic equiv- 
alent luminosity, while their origin still remains as a mystery. 
The redshift of one GRB 990123 was z = 1.6, and its total 
energy emitted as gamma-rays, estimated by its energy flux 
and redshift, is E ~ 3 x 10 54 erg assuming isotropic radiation 
(Kulkarni et al. 1999). This energy is equivalent to a rest 
mass energy of 1.7M Q c 2 and it is almost impossible to ex- 
plain by explosions of objects with stellar mass scale. This is 
why most researchers in this field now consider that GRBs are 
likely strongly collimated, with a typical collimation factor of 
47r/AS! > 100. It is then important to test this hypothesis by 
observations, for better understanding of the nature of GRBs. 

A direct consequence for such GRB collimation is that the 
rate of GRB afterglows may be much higher than those associ- 
ated with observable prompt GRBs. GRBs are believed to be 
produced by dissipation of kinetic energy of ultra-relativistic 
outflow from the central engine with a Lorentz factor of 
T ^100-1000. The outflow is eventually decelerated by in- 
teraction with interstellar matter, just like supernova remnants, 
to produce radiation called afterglows. After Y decreases to 
T - 0rj, where 6*j e t is the opening angle of the jet, the beam- 
ing of radiation is wider than that of outflow, thus the afterglow 
becomes observable from directions different from that of the 
original gamma-ray radiation. At the same time, the sideway 
expansion increases the opening angle of the jet from the origi- 
nal angle of 9j et . 



Therefore a serendipitous search of GRB afterglows without 
prompt GRBs is an important test for the beaming of GRBs 
(Rhoads 1997). Several efforts have been made in various 
wavebands. Past X-ray surveys set a constraint of b le j < several 
(Grindlay 1999; Greiner et al. 2000), where b m \ is the "rela- 
tive" beaming factor, i.e., the ratio of all afterglow rate includ- 
ing orphans to the rate of those with on-axis GRBs. 2 Perna & 
Loeb (1998) set a crude upper limit on b K \ < 1.5 x 10 3 by ra- 
dio source counts. Schaefer (2002) searched a 264 deg 2 field 
down to R ~ 21, but no afterglow candidate was found. The 
past survey for high-redshift supernovae should also give some 
constraint on the orphan GRB afterglow rate (e.g., Rees 1999). 
Recently, Vanden Berk et al. (2002) reported a possible or- 
phan afterglow found in the initial 1500 deg 2 data of the Sloan 
Digital Sky Survey (SDSS). However, Gal- Yam et al. (2002) 
reported that the host galaxy of the SDSS transient is a highly 
variable, radio-loud active galactic nuclei (AGN) and hence it 
was almost certainly not related to a GRB. 

The interpretation of these results is, however, not straight- 
forward. They must be compared with a realistic theoretical 
prediction of b re \ which is generally dependent on the waveband 
and sensitivity of surveys. Recently, Dalai et al. (2002) argued, 
based on a simple analytical investigation, that the value of b re \ 
is not sufficiently large for the sensitivity that can be achieved 
by the current facilities, even if GRBs are strongly collimated. 
They also argued that the rate of orphan afterglows is insen- 
sitive to the jet opening angle, and hence they concluded that 
the orphan afterglow search hardly constrains the GRB colli- 
mation in practice. What makes the situation even more com- 
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2 We used the term "relative" since b sc \ reflects the relative beaming of radiation between prompt GRBs and typical afterglows detectable in an orphan afterglow 
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plicated is the possibility that orphan afterglows can be pro- 
duced by other processes rather than the collimation; for exam- 
ple, a 'failed' GRB or a 'dirty fireball' may produce afterglows 
at longer wavelength while there is no prompt gamma-ray emis- 
sion (Huang, Dai, & Lu 2002). 

In this paper we present a realistic prediction of orphan af- 
terglow rate as a function of search sensitivity in various wave- 
bands, based on a popular afterglow model of collimated GRBs 
that has been tested against a number of afterglow observations. 
We find that the orphan afterglow rate rapidly increases with 
GRB jet collimation roughly as b K \ oc 6^ v if GRB jets have 
roughly constant energies against various 6*j et , as suggested by 
recent studies (Frail et al. 2001; Panaitescu & Kumar 2001). 
Therefore, the orphan afterglow search could still be powerful 
to get information on collimation of GRB jets. Even if the theo- 
retical prediction is also affected by possible other components 
of ejecta from successful or failed GRBs with lower Lorentz 
factor, it would be useful to make a prediction for the abun- 
dance of orphan afterglows by the 'nominal' (i.e., without any 
other hypothetical components) jet model of GRB afterglows, 
as realistic as possible, since the jet model is a relatively reli- 
able method of transforming the on-axis light-curves to off-axis 
ones, compared with other theoretical possibilities of orphan 
afterglows. It can be used as a baseline prediction when one 
interprets the results of past and future searches for faint extra- 
galactic transient objects. 

To make plausible predictions, we take the observed radio, 
optical, and X-ray emission of ten GRB afterglows with a 
wide variety of estimated jet opening angle — 970508, 9805 19, 
990123, 990510, 991208, 991216, 000301c, 000418, 000926, 
and 010222 — and calculate their emission as it would be seen 
at various angles relative to the jet axis. This is done within 
the framework of collimated jets decelerated by the circum- 
burst medium (Meszaros & Rees 1997) and undergoing lat- 
eral expansion (Rhoads 1999), using the treatment presented by 
Panaitescu & Kumar (2000) and Kumar & Panaitescu (2000). 
With the aid of this model, Panaitescu & Kumar (2002) have 
determined the physical parameters of the ten afterglows by 
modeling their broadband emission. We employ these ten sets 
of jet parameters to calculate their emission for any off-axis 
observer location. The detection rate as a function of the X- 
ray/optical/radio search sensitivity is further obtained by inte- 
grating over the viewing angle and over the GRB rate history in 
the universe. 

The paper will be organized as follows. §2 reviews the model 
of collimated afterglows used here, and some discussions are 
presented on the physics of collimated afterglows concerning 
the detectability of orphan afterglows. §3 is for the formulation 
of orphan afterglow rate calculation. The results are presented 
in §4; they are compared with various past search results and 
predictions for the current and future facilities are made. Dis- 
cussions are given in §6, including comparison with previous 
studies, another possible picture of GRB jets, and caveats of 
our predictions. Our main conclusions are presented in §7. 

2. THE AFTERGLOW MODEL OF COLLIMATED GRBS 

The dynamics of a GRB jet is calculated by tracking its total 
energy (allowing for radiative losses) and mass. The jet spreads 
laterally at the sound speed, and it is assumed to be uniform 
within its aperture (i.e. the same energy per solid angle in any 
direction) and to have sharp boundaries. The afterglow syn- 
chrotron emission and the radiative losses are calculated from 



the strength of the magnetic field in the shocked gas and the 
power-law distribution of shock-accelerated electrons, taking 
into account their radiative (synchrotron and inverse Compton) 
cooling. To obtain observer light-curves, the jet emission is in- 
tegrated over its surface, allowing for the differential relativistic 
beaming and the spread in the photon arrival time across the jet 
surface. 

The afterglow jet model has three parameters that determine 
the jet dynamics (the initial jet energy £j et , initial jet half-angle 
6*j e t, and external particle density n ext ), and three parameters per- 
taining to the microphysics of shocks (the fraction of the post- 
shock energy density in magnetic fields e B , the fractional energy 
of the least energetic injected electrons e e , and the power-law 
index p of the shock-accelerated electron energy distribution). 
For hard electron distributions with p < 2, two other parameters 
become relevant: the total fraction of the internal energy of the 
shocked fluid imparted to electrons, which sets a "cut-off" of 
the electron distribution at high electron energies, and the larger 
(than 2) power-law index of the electron index above the cut- 
off. The passage of the spectral break associated with the elec- 
tron distribution cut-off through the optical domain may yield 
a light-curve break, as suggested by the sharp break seen in 
000301c and the steep decay of the optical emission of 991208. 

Figures 1-3 show the X-ray, optical, and radio afterglow light 
curves for the above ten jets, as they would be seen from various 
viewing angles relative to the center of the jet, 6 b S . It can be 
seen that the behavior of off-axis light curves is very different 
for different sets of model parameters, even when the opening 
angle of the jet, 6>j et , is similar. As illustrated in Figure 2, 'the 
best case', i.e., the highest rate for orphan afterglow searches 
is expected for GRB 991216 (#j e t = 2.7°), whose peak of op- 
tical light curves is at R ~ 24 (when it is placed at z = 1) for 
an observer located with an angle 6> t, s = 30°. If one searches 
orphan afterglows with a sensitivity better than this, one could 
find them at more than 100 times higher rate than that expected 
from spherical GRBs (i.e., b m \ > 100). The other extreme is the 
cases of GRB 980519 or 990123. Even though their jet open- 
ing angle #j et is not much different from that of GRB 991216, 
the peaks of optical light curves are difficult to detect when an 
observer is located at o b s > 10°. In these cases the relative 
beaming factor is at most b m \ < 10. 

The dependence of b K \ on the jet initial opening (and other jet 
properties), as well as the reason for which b K \ is maximal for 
the 991216 afterglow, can be obtained using a simplified model 
where all the afterglow emission originates from a single point 
on the jet axis (Dalai et al. 2002, Granot et al. 2002) whose dy- 
namics is that of an expanding jet (Rhoads 1999). For narrow 
(#jet -C 1 rad), relativistic jets and a homogeneous medium, it 
can be shown that the jet Lorentz factor 7 is 

7( ° } "M^J { n= 1/2, t >t Jfl ' (1) 
where the subscript "0" indicates the time (since the prompt 
gamma-ray emission) measured by an on-axis observer (9 b S = 
0), tjfi being the "jet-break time", i.e. the time when 7 = 8~} t 
and the on-axis observer "sees" the entire jet surface and a 
break in the light curve. Taking into account that the photon 
arrival time for an arbitrary observer location 6> b s is given by 
dt/dto = 2j 2 (l -/3cos# bs), fi being the jet speed in units of 
the speed of light, equation (1) gives the jet Lorentz factor as a 
function of the time t when photons arrive at o b s : 

7(0 ^el / t 



/"(7) f»(fi£) W^obs) 



(2) 
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where 

/( 7 ) = cos0 obs + 8[2-ln(min{l, 7 jet })] 7 2 sin 2 (0 obs /2) . (3) 

The above /( 7 ) relates the photon arrival time at 9 ^ to that 
for an on-axis observer: t = /( 7 )fo, therefore, in equation 
(2), tj(8 ohs ) = f(0~} t )tj,o. For an observer located well out- 
side the jet opening (6 \, s 3> 6*j et ) the received jet emission is 
Doppler boosted in frequency by the same factor T> = [ 7 (1 - 
/3cos6* bs)] _1 , its intensity being relativistically enhanced by 
T> 3 . Thus the received afterglow flux F(v : t) at # bs is related 
to that seen by an on-axis observer, Fo(*Vo)> through 

F{u,t) = g-\ 1 )F a {gv,t/f), (4) 

where 

g( 7 ) = cos6» obs + 4 7 2 sin 2 (6l obs /2) . (5) 

The observed GRB afterglows have power-law optical 
spectra and light-curves: Fo(zVo) cx v~ 8 t^ a (for the above- 
mentioned ten afterglows, S G [0.6, 1.5], a = ct\ £ [0.7, 1.7] at 
to < 0,0 > and aEa 2 £ [1 .6, 3.0] at to > f/,o). Therefore equation 
(4) leads to 

F(^,f)cxr( 7 )[ g ( 7 )]- (3+5) r Q , (6) 

with 7 (r) given by equation (2). 

For t < tj, when 7 0j et > 1, the above equations yield 7 cx r 3 / 2 



and F(v,t) cx 7 



2(ai-(5-3)+-a 



r a, ^ f 3(3+5)-4 Ql in the 6 ) obs ^ : j imit ^ 



i.e. a light-curve with a sharp rise. At t > tj and before the time 
t p when 7 (f /) ) = 9~l s , 7 oc exp{-2f/f/} and F{v,t) cx exp{4(3 + 
(5- a2)(///j)}^ thus the light-curve continues to rise. At t > t p 
the / and g functions asymptotically approach unity, so that 
the evolution of the jet Lorentz factor and the afterglow light- 
curves become those for an on-axis observer: 7 cx f" 1//2 and 
F(v : t) cx t~ ai , respectively. Therefore the afterglow light-curve 
for 6* bs peaks around the time t p when j(t p )9 bs = 1 • Taking into 
account that t p = f[j(tp)]tp.o, where t pfl ~ (0 O bs/0jet) 2 O,o (from 
eq. [1]) is the corresponding photon arrival time for an on-axis 
observer and f(9~ bs ) ~ 5 + 21n(# b S /#j e t) in the 6* b s -C 1 limit, 
it follows that the afterglow light-curve seen by an observer at 
O bs peaks at 

2 

0,o • (7) 



t p = [ 5 + 21n- 



i/jet 



In the same 9 obs <C 1 limit, g[ 7 (f p )] : 
and (7) give for the peak flux 



2, therefore equations (4) 
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It may not be easy to infer t p for an observed orphan afterglows, 
since we do not know the precise time of the prompt burst time. 
However, if there are enough number of data points to construct 
a light curve, it may be possible to infer t p only from orphan af- 
terglows. Then, together with an empirical relation between the 
jet-break time tj$ and flux Fo(v,tj t o), calibrated with the after- 
glows which were seen on-axis (i.e. preceded by a GRB), equa- 
tions (7) and (8) can be used to determine the relative observer 
location O bs/#jet> tj$, and Fq{v, tj.o) from the inferred peak time 
t p and observed flux F{v,t p ) of an orphan GRB afterglow. 

From equation (8) it follows that only jets seen at an angle 
less than 

l/(2a 2 ) 



o-p+ ^o^O.o) 

F\i m 
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jet •> 



(9) 



can be detected above a given detection threshold Fn m . 
Therefore the relative beaming factor is b K \ = (9 m - dx /9j et ) 2 cx 



[2 _<5 Fo(^,0.o)] 1,/Q2 - The b K \ that can be inferred for each af- 
terglow shown in Figure 2, for a given threshold, is consistent 
with this estimation of b m \. Among the set of ten afterglows, 
991216 has one of the largest jet-break fluxes F a (i>,tjxi) (note 
that for this afterglow f,.o < 1 day, the sharper break seen in 
Figure 2 at several days being due to the passage through the 
optical band of a spectral break), hardest optical spectra, and 
shallowest decays after the jet-break time. For these reasons, its 
b m \ is the largest in our sample of afterglows, and GRB991216 
dominates the detection probability of orphan afterglows. 



3. FORMULATIONS FOR THE RATE CALCULATION 

First we calculate the expected number of orphan afterglows 
per unit solid angle N ex? (Fii m ), which can be detected by a snap- 
shot observation with a given sensitivity Fn m . This quantity di- 
rectly provides the detectability of afterglows in a search where 
transient objects are found by comparison of several snapshot 
images taken at time intervals longer than the duration T over 
which the afterglow is brighter than the sensitivity. On the 
other hand, the detection rate in a survey of consecutive mon- 
itoring with a time scale longer than T can be estimated by 
^ex P ~ N ex? /T per unit solid angle and observation time. 

The afterglow flux at any frequency and observer's time 
viewed by an observer with arbitrary 9 Q \, S is calculated with the 
model presented in the previous section. From there we can 
calculate the time duration for an observer, T(x,z, O bs,^iim), 
during which an afterglow is brighter than a given sensitivity, 
where z is redshift and x symbolically represents a specific type 
of GRB afterglows. Then, the expected number N exp per all sky 
can be written as: 



N exp (Fi im ) = JdxJdzJ d9 ohs 

dV RcRB(x, z) 
dz l+z 



sniffs 



-T(x,z,0 O b S ,Fii m ) 



(10) 



where (dV /dz) is the standard comoving volume element in all 
sky, and /?grb(^,z) is the comoving GRB rate density of type x. 
It should be noted that /?grb is the true rate of GRBs including 
unobservable ones whose jet are not directed to us. 

As shown in the previous section, the detectability of orphan 
afterglows sensitively depends on afterglow parameters. In or- 
der to take this into account and make the most plausible pre- 
diction in an empirical way, we replace the integration over x 
by the sum of the ten well-observed afterglows presented in 
the previous section. We assume that all the population of ob- 
servable GRBs is represented by those ten bursts with an equal 
weight, and it is independent of redshift. If detection of GRBs 
is flux limited, then a simple sum of the ten observed bursts hav- 
ing very different gamma-ray luminosities might induce some 
bias from the true population. However, GRBs are generally 
bright enough to be detected even at very large cosmological 
distances, and hence the effect of flux limit is expected to be 
insignificant. In addition, the redshifts of the ten GRBs are 
in a rather small range. Therefore, we expect that the above 
treatment is not unreasonable. On the other hand, it should be 
noted that we must take into account the estimated opening an- 
gle of jets; the GRB rate inferred from the observed number of 
GRBs assuming isotropic emission (R lm ) is not the true GRB 
rate, but it should be multiplied by the inverse of sky coverage 
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of gamma-ray emission, 2/(1- cos 6*j et ), to obtain /?grb ■ 3 Then, 
the integration over x in equation (10) should be replaced by the 
sum of JVgrb (=10) GRBs, as: 



dV R iso (z) 
x— , 2 



/ l-cos0 jet A 1 

I ^ J UiZj^obsj^lim) • 



(ID 



It should be noted that the ten GRBs chosen by Panaitescu & 
Kumar (2002) are well observed, with good temporal coverage 
at various frequencies, and showing a light curve break. This 
suggests that there may be a systematic bias toward strongly 
collimated GRBs in our selection. We note that the ten GRB 
afterglows used here represent about a third of all bursts for 
which optical afterglows were found. This indicates that the 
formulation above overestimates by a factor of up to three the 
real orphan afterglow rate if the other 2/3 GRB afterglows are 
isotropic and hardly contribute to N exp . This is probably not 
much larger than other model uncertainties. Therefore we be- 
lieve that the ten GRBs represent the whole GRB population 
reasonably well. 

We assume that Ri so (z) traces the cosmic star formation his- 
tory (Totani 1997; Wijers et al. 1998). Based on the recent 
observational estimates (see, e.g., Totani & Takeuchi (2002) 
for a summary of recent studies on this issue), we assume that 
Riso(z) oc (1 +z) 3 - 77 at < z < 1, and R iso (z) is constant at z > 1, 
when the Einstein-de Sitter universe is adopted. We set an up- 
per cut-off for redshift as z u = 5. The normalization of this rate 
density is fixed at z = 0, as R iso (0) = 2 x I0~ w (h /0.7) 3 yr-'Mpc" 3 
where h = Hq/( 100km s^Mpc" 1 ), based on the fitting to the ob- 
served flux distribution (i.e., log/V-logP relation) of the BATSE 
GRBs (Totani 1999; Schmidt 1999). Throughout this paper we 
adopt a standard A-dominated universe with (h, firj, Qa) = (0.7, 
0.3, 0.7) and the cosmology-dependence of /?i S0 (z) is corrected 
appropriately. 

4. RESULTS 

Figs. 4-6 show the result of calculation for 7V exp , as a function 
of sensitivity in X (1 keV), optical (R), and radio (5 GHz) wave- 
bands, respectively. We also show the rate expected only from 
on-axis afterglows, i.e., associated with prompt GRBs. This 
rate, N on , is obtained by replacing the range of integration over 
#obs = - 7r by 6> bs = — 0jet in eq (11). The relative beaming 
factor, b m \ = N exp /N on is shown in the lower panel. The mean 
values of redshift and T are also shown; here we calculated av- 
erage of logz and log T because of considerable scatter of these 
values at a fixed sensitivity. We also give another mean time, 

1_\ = ([dN^/£ 

J / \ N exp 

where integration with dN exp represents that of eq. (11). We 
give this quantity because we can translate N exp into detec- 
tion rate per time in a consecutive monitoring observation by 
R exp = N ex? (l /T), taking into account the dispersion of T. For 
comparison, N exp for type la and II supernovae calculated by 
Woods & Loeb (1998) for the optical band is also shown in Fig. 
5. 

A clear trend can be seen in all wavebands; the relative beam- 
ing factor b m \ is rapidly increasing with the sensitivity limit, 

3 Here, we have assumed that the opening angle of gamma-ray emission is the same with that of the jet inferred from the afterglow fitting, but this is not necessarily 
true. We will discuss this point later in §6.3. 
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while the mean redshift is not. These results reflect the fact that 
intrinsically faint afterglows at small distances with large 9 b s 
become dominant when the search sensitivity is improved. For 
orphan afterglow search with b K \ > 10, the required sensitivi- 
ties are vF v ~ 10~ 14 -10~ 13 erg cm~ 2 s _1 and R ~ 22, in X-ray 
(IkeV) and optical bands, respectively. In the radio bands, b re \ 
exceeds 10 in all flux range shown here, while a sensitivity of 
~ lO^Jy is necessary to increase b K \ to more than 100. 

Figure 7 shows the contribution from each of the ten GRBs to 
the total rate of orphan afterglows. As expected, GRB 991216 
dominates in most cases, while others are important as well in 
bright flux range in radio bands. It is interesting to check how 
the prediction is changed when the dominant GRB 991216 is 
removed in the sample for averaging. In the optical band, the 
mean rate is reduced by a factor of 2.2 and 4.5 at the search 
sensitivities of R = 20 and 25, respectively, when GRB 991216 
is removed. 

5. COMPARISON WITH PAST SEARCHES, AND PROSPECTS 
FOR FUTURE PROJECTS 

In this section we calculate the number of orphan afterglows 
expected in several past surveys and compare them with the 
reported results. We also make predictions for the number of 
orphan afterglows expected in future surveys. The summary of 
our results is given in Table 1 . 

5.1. X-ray Observations 

There are two papers which constrained the orphan after- 
glow rate in X-ray band: Grindlay (1999) by the Ariel 5 sur- 
vey with a sensitivity of ~ 10~ 10 ergs cm~ 2 s _1 (2-10keV), and 
Greiner et al. (2000) by the ROSAT All Sky Survey (RASS) 
which is sensitive to ~ 10" 12 ergs cm" 2 s _1 (0.1-2.4keV). Fig- 
ure 4 shows that the relative beaming factor b K \ is about 1 .6 and 
3.0 for the sensitivities of the Ariel 5 and RASS, respectively. 
Considering this small number, the negative result of Grindlay 
(1999) seems consistent with our expectation. The exposure of 
RASS is 76,435 deg 2 days, and using (1 /T)~ l ~ 0.8 day at the 
RASS sensitivity, this is equivalent to a snapshot observation of 
~ 94, 000 deg 2 . Our model predicts about 8 GRB afterglows for 
this sky coverage, among which ~3 are expected to be on-axis. 
Greiner et al. found 23 candidates of afterglows, but they ar- 
gued that the bulk of these are likely to be nearby flaring stars. 
If one removes those which are suspected as flaring stars, the 
remaining candidates are at most ~10. It is interesting that this 
number is very close to our expectation. 

The advanced satellites such as Chandra or XMM-Newton 
have a typical sensitivity limit of ~ 10~ 15 erg cm" 2 sec _1 for 
point sources, with a field of view of ~ 10 3 arcmin 2 . At such 
a sensitivity, the relative beaming factor is increased to b m \ ~ 
40. Our model predicts that the probability of finding an or- 
phan afterglow in the field of view of one snapshot observation 
is ^0.02. This is apparently small, but accumulation of archive 
data might be useful to search and constrain orphan afterglows 
in the future, though discrimination from flaring stars could be 
a major problem again. 



5.2. Optical Observations 
5.2. 1 . Past Optical Surveys 
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Schaefer (2002) searched orphan afterglows in a 264 square 
degree field with a sensitivity of R = 21. The field was exam- 
ined nightly, with a total duration of 33 days. He found no af- 
terglow candidates. By using (1 /T)~ l ~ 5.1 days of our predic- 
tion, this observation corresponds to a snapshot observation of 
264 x 33/5.1 ~ 1,500 square degree. Our model predicts that 
about two orphan afterglows are expected, which is b K \ ~ 6.2 
times higher than the case of no beaming. Therefore, our model 
is marginally consistent with his result. On the other hand, 
there is a good chance to detect orphan afterglows by contin- 
ued search and the effort in this direction is encouraged. 

High redshift supernovae have been intensively searched 
down to R ~ 23, for the purpose of determination of cosmo- 
logical parameters (e.g., Schmidt et al. 1998; Perlmutter et 
al. 1999). The total exposure is about a few tens of "square 
degree years" (Rees 1999). They are normally searching su- 
pernovae with two images separated by one month, which is 
longer than the characteristic afterglow time scale at this sensi- 
tivity ((1/r)" 1 <~ 18days). Therefore, the exposure is equiva- 
lent to 12 times searches in a field of a few tens of square de- 
gree. Then we estimate the expected detection number is about 
1.6 in the past, which is b K \ ~ 14 times higher than the case of 
no beaming. It should be noted that the past supernova searches 
(including that of Schaefer) are not perfectly suitable to orphan 
afterglow search. Supernovae are found only in one bandpass 
(i.e., no color information), and spectroscopic follow-up is done 
only for a small part of them showing good properties for the 
cosmological use. It should also be noted that they sometimes 
find "mystery objects", which decay faster than supernovae and 
have no host galaxies (Schmidt et al. 1998). These objects 
might be orphan afterglows. 

5.2.2. Sloan Digital Sky Survey 

Recently Vanden Berk et al. (2002) reported an interesting, 
highly luminous transient object which might be an orphan af- 
terglow, found in the initial 1500 deg 2 data of the Sloan Digital 
Sky Survey. However, Gal- Yam et al. (2002) claimed that the 
host galaxy of the SDSS transient is in fact an unusual radio- 
loud AGN showing strong variability. Therefore, it seems very 
likely that the SDSS transient was an AGN flare, though the 
possibility of an orphan afterglow cannot be completely re- 
jected. Here we give an estimate of expected number of or- 
phan afterglows in the SDSS data taken so far. Vanden Berk 
et al. (2002) searched transient objects by the flux difference 
between imaging observations and later spectroscopic observa- 
tions. Although the SDSS imaging observation has a sensitivity 
of r' ~ 23, spectroscopy is not performed for all objects; the ob- 
ject was selected as a quasar candidate by its colors. For low 
redshift quasars whose colors are similar to GRB afterglows, 
the sensitivity is ~ 19 (Stoughton et al. 2002). Since R-i' - 
for typical GRB afterglows, we can consider that the sensitivity 
of the search made by Vanden Berk et al. is R <~ 19, with a sur- 
vey area of 1500 deg 2 . Then, we found that the expected num- 
ber of orphan afterglows to be <~ 0.2, which is b K \ ^3.0 times 
higher than the case of no beaming. This number is small, but 
not extremely small compared with one detection. (Note, how- 
ever, that the number becomes smaller by a factor of 10 when 
we apply the brightness of the SDSS transient, R ~ 17, as the 
search sensitivity.) Our model predicts the mean redshift at this 
sensitivity to be 10^ log ^ ~ 0.6, and the la dispersion around 
this mean is <ri ogz ^0.5. Therefore the redshift z = 0.385 of the 
host galaxy well falls in the plausible range. The characteris- 



tic time scale predicted by the model, lO^ 7 ") ^3.5 days with 
ciogr ~ 0.7 is also consistent with the modest variability shown 
by the SDSS objects during the first two observations separated 
by 2 days. Therefore, it is not unreasonable even if the SDSS 
transient is an orphan afterglow, although the strong variabil- 
ity of the host galaxy indicates that the transient is much more 
likely an AGN flare. More statistics is obviously needed for 
stronger conclusions. 

When the SDSS project is completed, the covering area in 
the northern sky will be increased to 10,000 deg 2 . Then we 
expect about 1.3 orphan afterglows, which seems still small to 
test our model with sufficient statistics by the future data. On 
the other hand, the SDSS has another observing mode in the 
southern sky, observing a 225 deg 2 field repeatedly to achieve a 
much deeper sensitivity limit than the northern sky (e.g., Ivezic 
et al. 2000). These southern data can be used for transient ob- 
ject search with a sensitivity of R <~ 23, which is significantly 
deeper than the search using spectroscopic data. The southern 
sky field is typically observed about four times in a year with 
a separation longer than one month. Then, at the end of the 
project after the planned five year operation, we expect an ef- 
fective survey area of 225 x4x5~ 4500 deg 2 (Z. Ivezic 2002, 
a private communication). Then we expect about 36 orphan 
afterglows, which is b re \ ~ 14 times higher than the case of 
no beaming. This number suggests that we may detect a sta- 
tistically meaningful number of orphan afterglows. Although 
spectroscopic information will not be immediately available, 
follow-up observations for host galaxies will give redshift in- 
formation. The five-bands photometry data will also help to 
discriminate the orphan candidates from other transient objects 
(Rhoads 2001). Contamination of highly variable AGNs, as in 
the case of the SDSS transient of Vanden Berk et al. (2002), 
can be removed by close examination of host galaxies in X- 
ray and/or radio bands and past records. (It should be noted 
that the host galaxy of the SDSS transient reported by Vanden 
Berk et al. has a strong radio emission that cannot be explained 
by star formation activity.) Finally, a significant part (~ 10- 
40%) of the northern sky will also be observed more than two 
times, because of overlaps of field-of- views. The time intervals 
of these repetitions are not simple but depending on the survey 
schedule. These data can also potentially be used for orphan 
afterglow searches. 

5.2.3. OGLE III 

The Optical Gravitational Lensing Experiment in the third 
phase (OGLE III, Udalski et al. 2002) has started its operation, 
and during the LMC and SMC season, it will cover a 85 square 
degree field every night or every second night for half a year. 
The limiting magnitude is about /= 19.5 at S/N=10 for a typ- 
ical exposure of 2 minutes, corresponding R ~ 20 for typical 
afterglow spectrum. Using (1/T)~ l = 3.0days, the total expo- 
sure of half a year is equivalent to a snapshot observation of 
85 x 180/3.0 ~ 5000 square degree. Then, about 2 or 3 orphan 
afterglows are expected. Since the variability time scale is only 
a few days at this magnitude, continuous monitoring of OGLE 
III with the time interval of one or two days is very useful for 
an afterglow search. The mean redshift is \0 l ° sz = 0.65 with a 
dispersion of <J\ ogz ~ 0.5. This means that a significant fraction 
of orphan afterglows detectable by OGLE III should have large 
redshifts of z > 1. Then the brightness of R ~ 20 would be 
much brighter than any kind of supernovae and hence orphan 
afterglows and supernovae can be discriminated. 
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5.2.4. ROTSE-III 

The ROTSE group is considering to use their ROTSE-III tele- 
scope for an orphan afterglow search, whose field of view is 
3.5 deg 2 (Smith et al. 2002; Kehoe et al. 2002). A search by 
planned four instruments will cover about 1400 deg 2 with a lim- 
iting magnitude better than 19th each night (C. Akerlof 2002, a 
private communication). This means that ROTSE-III has a ca- 
pability of doing an orphan afterglow search equivalent to that 
made by SDSS (Vanden Berk et al. 2002) only in one night. 
If such a search is continued for a year, the effective sky cov- 
erage could reach ~ 1400 x 365/1.3 ~ 3.9 x 10 5 deg 2 , taking 
into account that (1/T)~ l = 1.3 days at this sensitivity. Then 
we expect more than 50 orphan afterglows with b re \ ~ 3. Since 
T < 1 day, more than one observations in a night would also be 
favored for sufficient time resolution. 

5.2.5. Subaru/Suprime-Cam 

The Subaru Prime Focas Camera (Suprime-Cam) of the 8.2m 
Subaru telescope is a unique facility with a field of view (FOV) 
of 30' x 30', which is more than 100 times wider the than typi- 
cal field of views of 8m-class telescopes or HST Therefore this 
instrument is the most suitable for an orphan afterglow search 
at the deepest sensitivities we can achieve. A sensitivity limit 
of R ~ 26 for point sources 4 is achieved by about 10 minutes 
exposures, and then it is possible to observe about 10 FOVs 
in a night with multi-band photometry. A clear advantage of 
the deep sensitivity is that we expect a large relative beaming 
factor, b K i ~ 50. Then we expect about 0.4 orphan afterglows. 
Therefore a detection is not extremely difficult when system- 
atic searches are performed using several nights. We show the 
expected number assuming a search over a 5-deg 2 field in Ta- 
ble 1 . If detected, a large value of b le \ strongly argues for the 
existence of many more orphan afterglows than those associ- 
ated with GRBs. At this sensitivity, the time scale of afterglow 
variability is increased to 10 logr ~ 150 days with ai og T ~ 0.4. 
Therefore, a longer time scale than supernova searches may be 
favored. 

A problem is discrimination from other transient objects. A 
few supernovae are typically found in one FOV of the Suprime- 
Cam down to / ~ 25 (M. Doi & N. Yasuda 2002, private com- 
munications). At this magnitude, orphan afterglows may not 
be sufficiently brighter than supernovae, while the SDSS tran- 
sient of Vanden Berk et al. (2002) was about 100 times brighter 
than the brightest supernovae. Most supernovae have thermal 
spectra that are curved compared with that of power-law GRB 
afterglows. Therefore, color-color plots by multi-band photom- 
etry can be used for discrimination (Rhoads 2001), though it is 
still uncertain whether this method can remove all supernovae 
including peculiar ones. Offset of optical transients from cen- 
ters of galaxies, close examination of host galaxies, and past 
records will be useful to remove AGNs, as mentioned in the 
previous section. 

5.2.6. GAIA 

The astrometric satellite GAIA will survey all the sky many 
times with a sensitivity of R ~ 20. Each location of the sky 
will be observed about 40 times separated by more than one 
month, during the whole mission (L. Eyer 2002, private com- 
munication). Therefore we expect about 720 orphan afterglows. 

4 This estimate is not seriously affected even if an afterglow is found in a host g: 
observation is limited by the sky background, and in most cases the surface brigh 

5 http://www.dmtelescope.org 



While the relative beaming factor is not large (b K \ ~ 4.3 at this 
sensitivity limit), the enormous number expected might be use- 
ful for statistical analysis of orphan afterglow rate. Because of 
small b le i, the cross check between gamma-ray observation is 
crucial to discriminate orphans from ordinary afterglows asso- 
ciated with observable GRBs. It is highly desired that a GRB 
satellite covering a significant part of all sky is working at the 
time of GAIA project. 

5.2.7. Dark Matter Telescope 

Even more powerful searches than that by the 
Subaru/Suprime-Cam would become possible at the sensitivity 
of 8m-class telescopes, by the planned Dark Matter Telescope 
(DMT) 5 , having a 7-deg 2 field of view. In the planned all-sky 
survey mode, the DMT will cover 20,000 deg 2 down to 24th 
magnitude twice in the same month, in which about 500 orphan 
afterglows are expected with b K \ ~ 20. On the other hand, the 
deep probe mode will cover ten 100-deg 2 fields with a sensitiv- 
ity of 29th magnitude. Then, in principle, a few thousands of 
orphan afterglows could be detected with b m \ ~ 200, although it 
is not yet clear how they can be discriminated from other tran- 
sient objects at such a deep sensitivity level. It should also be 
noted that the characteristic time scale T is more than 500 days 
for this sensitivity, and hence a sufficiently long time interval is 
required for the search. 

5.3. Radio Observations 

Perna & Loeb (1998) used faint source counts in radio bands 
to constrain the orphan afterglow rate. The radio source counts 
at 8.44GHz is - 3 x 10 6 in all sky for S > 0.1 mJy (Windhorst 
et al. 1993; Becker, White, & Helfand 1995). Following Perna 
& Loeb, we assume that about 3% of sources at this sensitivity 
are variable, yielding an upper bound on the number of orphan 
afterglows in all sky as < 9 x 10 4 . The power index of the radio 
spectrum is changing from the initial value of S = -1/3 (be- 
low the peak frequency) to S ~ 1 (above the peak frequency), 
but typically it is S ~ 0.5 at a few hundreds days after the 
burst. Taking into account this correction from 8.44 to 5GHz, 
we obtain ~2x 10 3 orphans expected by our model. There- 
fore the current upper bound on variable radio source counts 
hardly constrains our model. In fact, the relative beaming fac- 
tor expected at this sensitivity is b K \ ~ 15, which is not greater 
than the maximum expected by the deepest optical search by 
Subaru/Suprime-Cam. One might have thought that radio band 
is potentially the best waveband for orphan afterglow searches, 
since radio afterglows are visible long after the bursts and hence 
we expect a large b m \. However, in order to achieve b m \ > 100, 
the sensitivity must be better than ~- 1—10 ^Jy. 

The difficulty in radio searches is discrimination from other 
variable sources such as AGNs. It is difficult to do this discrimi- 
nation only in the radio band, but might be possible if we utilize 
information in other wavebands. Intensive examination of vari- 
able objects found in the radio band with cross-checks between 
the radio, optical, and X-ray surveys may allow us to find an or- 
phan afterglows in a large (~ 40-50) sample of variable AGNs 
at sensitivities of ~ 0.1 mJy. 

A unique characteristic of radio afterglow, which is different 
from X-ray or optical bands, is that b K \ does not decrease but 

axy whose luminosity is brighter than the afterglow, since generally ground-based 
less of galaxies is not brighter than the sky. 
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stay roughly constant at b K \ > 10-20 when the search sensitiv- 
ity is decreased to F„ > lmJy (brightest flux region). This is 
coming from the properties of the radio light curves; the peak 
of radio light curves occurs at relatively late time even in the 
on-axis case. Therefore, the radio afterglows cannot be de- 
tected at large cosmological distances even when an observer 
is located on the jet axis, while on-axis early afterglows in X- 
ray or optical bands are very bright and they can be detected at 
almost all the cosmological distance scale. As can be seen in 
Fig. 3, the peak radio flux is almost independent of o \, s in some 
dominant GRBs such as 991216, 000301c, and 000926. There- 
fore b K i converges into a finite value (~10-20) with Fn m — > oo. 
On the other hand, as mentioned earlier, very bright on-axis 
afterglows at large cosmological distances are dominant in the 
brightest flux range in the X-ray or optical bands, making b m \ 
to converge to unity with F\ lm — > oo. The sharp decrease of 
the mean redshift (z) of radio afterglow with increasing flux in 
a flux range of F v > lmJy, where b K \ is roughly constant, is 
also consistent with this interpretation (see Fig 6). This effect 
then suggests that a relatively shallower, but wide-field search 
is an efficient way to constrain the orphan afterglow rate in ra- 
dio bands, rather than deeper and narrower searches. In such a 
survey, low-redshift afterglows with large T and o \, s should be 
dominant, which may be called as GRB remnants rather than 
afterglows (Paczyhski 2001; Ayal & Piran 2001). 

Levinson et al. (2002) made a search of orphan afterglows by 
comparing the FIRST and NVSS surveys. Their sky coverage 
is 5990 deg 2 , and they found 26 candidates of orphan after- 
glows, for which they argued that these are unlikely to be ra- 
dio supernovae, while the possibility of radio-loud AGNs can- 
not be rejected. The number in 5990 deg 2 may be increased 
to ~ 65 when corrections for incompleteness are made. Their 
search sensitivity is 6 mJy at 1.5 GHz, and hence we trans- 
formed this into 3.3 mJy at 5 GHz band again assuming 5 = 0.5. 
Then we found that our expectation for this search is 7Y exp ~ 2.0 
at b K \ ~ 16. Further observational inspection of these candi- 
dates is necessary, and such effort may reveal some orphan af- 
terglows. If, instead, the majority of the candidates turned out 
to be orphans, it would indicate a very high orphan afterglow 
rate, which cannot be explained within the jet model with sharp 
edges. 

A future project, Allen Telescope Array (ATA), for the SETI 
has a sensitivity of about 0.3mJy at 1.4GHz (0.16mJy at 5GHz 
if 8 = 0.5) in one minute of integration time, with a pixel size of 
~ 1 arcmin diameter and a total beam area of 2.5 degree diam- 
eter. At this integration time, it could cover 25% of the north- 
ern sky everyday (L. Blitz 2002, a private communication, see 
also http://www.seti-inst.edu/science/ata.html). Although the 
sky coverage is not much better than the FIRST/NVSS sur- 
veys utilized by Levinson et al. (2002), the better sensitivity 
increases the expected number of orphan afterglows greatly to 
~ 200, out to the redshift of z ~ 0.6 with a characteristic vari- 
ability time scale of 200-300 days. Furthermore, ATA will pro- 
vide radio light curves sampled everyday, which would be very 
useful to check whether a transient source is an orphan after- 
glow or not. 

6. DISCUSSION 

6.1 . Comparison with Previous Work 

Our result that the detection rate of orphan afterglows of col- 
limated GRBs in optical bands could be much higher {b K \ > 100 
at R ~ 24 for GRB 991216) than in the case of spherical GRBs 



seems apparently in contrast to that of Dalai et al. (2002), who 
showed a result that the detection rate of orphan afterglows is 
insensitive to the jet opening angle 0j et , and b m \ is constant at ~ 
3-4 for a search with a sensitivity of R ~ 27. We note that the 
result of Dalai et al. (2002) was derived assuming fixed values 
for the jet-break time and luminosity at the break time for an 
on-axis observer. In this case, the maximum viewing angle at 
which the afterglow can be detected [0 max in eq. (9)] is propor- 
tional to #j et , and hence the relative beaming factor b re \ does not 
depend on #j et , as can be seen from equation (9). 

However, if the energy of afterglow jets does not vary much 
among GRBs, as suggested by Frail et al (2001), despite a wide 
range of initial jet angles 9 iet (Panaitescu & Kumar 2001), the 
flux at the jet-break time Fo(v,tjfi) measured by an on-axis ob- 
server is strongly dependent on 6- ]eX , thus, far from the assump- 
tion of constancy made by Dalai et al (2002). Within the frame- 
work of relativistic jets, it can be shown that at optical frequen- 
cies Fo(v,tjfi) oc 9J^ P , where p is the index of the power-law 
electron distribution. Furthermore, after the jet-break time the 
slope -«2 of the afterglow decay is —p (Rhoads 1999). Then 
equation (8) leads to a peak flux F(v,t p ) for an off-axis ob- 
server that is independent of 6- ]eX . This conclusion can also be 
reached by noting that the jet dynamics after the jet-break time 
is independent of 0j et (Granot et al. 2002). As shown in §2, for 
an off-axis observer the afterglow light-curve peaks after the 
jet-break time, therefore the peak flux, which is determined by 
the jet dynamics and observer location, does not depend on the 
initial jet opening. That F(u,t p ) is independent of 9j et implies 
that 9 m - dx does not depend on #j et and the original expectation 
&rei oc (0 m . dX /9j et ) 2 oc 6>j" e 2 is restored. 

The above result regarding the constancy of F(v,t p ) is illus- 
trated in in Figure 8, where we show the light curves of optical 
afterglows with various values of jet for an observer located at 
#ob s = 20°. Here we used a fixed set of typical afterglow model 
parameters (other than 8j et ): £j et = 3 x 10 50 erg, n ext = 1cm" 3 , 
e e = 0.05, e B = 3 x 10" 3 , and p = 2. The peak flux hardly changes 
with #j et at #j et < 5°. Figure 9 shows the beaming factor b le \ as a 
function of #j et for various sensitivities. Note that b K \ becomes 
larger for more collimated GRBs, asymptotically reaching the 
expected relation b m \ oc 9^ t when 9 jet — > 0. Thus, orphan after- 
glow searches could give useful information on the GRB colli- 
mation. 

6.2. Uniform versus Universal Jet Profile 

In this work we have assumed a conical jet with a sharp edge, 
a uniform energy per solid angle within the jet opening, and no 
energy outside (the uniform jet model). Such a model is ap- 
propriate if the angular distribution of the energy has a char- 
acteristic angular scale and decreases rapidly beyond it, e.g., 
an exponential profile oc exp(-9/8j et ). In this model, the ob- 
served anticorrelation between isotropic equivalent luminosity 
and jet break time can arise if jets have roughly the same en- 
ergy but different opening angles 0j et among bursts (Frail et al. 
2001). However, if the energy per solid angle has a large vari- 
ation but no characteristic scale (e.g. a power-law distribution) 
a completely different picture is possible. Some recent papers 
proposed that the observed luminosity - break-time anticorre- 
lation can be explained by a universal jet with a non-uniform 
profile, observed at different viewing angles (Rossi et al. 2002, 
Salmonson & Galama 2002, Zhang & Meszaros 2002), which 
can be an alternative to the uniform jet model. In such a model, 
the angular distribution of jet energy per unit solid angle should 
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be a power-law (d£/dCl) oc 9~ 2 to reproduce the above men- 
tioned correlation. 

It is not straightforward to predict how the orphan afterglow 
rate is changed when such a picture is adopted rather than the 
uniform jet model. If we have an ideal gamma-ray detector that 
can detect all GRBs everywhere in the universe, then we ex- 
pect that b K i cannot be much greater than unity in the universal 
jet picture, because of the following reason. We expect orphan 
afterglows only when viewing angles larger than an angle f? 7 
corresponding to either the maximum angular spread of the jet 
or the angle at which small energy per solid angle and/or the 
Lorentz factor yield a barely detectable gamma-ray emission. 
The orphan afterglow rate should then be close to that predicted 
by our model for a jet of opening angle 6 1 . An estimate of 9 1 
can be obtained as following. The isotropic equivalent gamma- 
ray energies calculated by Bloom, Frail & Sari (2001) for bursts 
with known redshifts span 2-3 orders of magnitude. If the GRB 
output is mainly determined by the jet energy per solid angle 
toward the observer, then the above distribution (d£ /dVl) oc 9~ 2 
implies that the dimmest GRB jets are seen at an angle at least 
10 times larger than that for the brightest jets. The latter angle 
should be around the smallest jet opening angle of 2° found by 
Panaitescu & Kumar (2002) by modeling the broadband emis- 
sion of ten afterglows using the uniform jet model 6 , leading 
to 9 1 > 20°. Since the relative beaming factor decreases with 
the jet opening angle, as shown in Figure 9, we expect that for 
the universal structured jet model b re \ <~ a few for the reason- 
able search sensitivities of R < 27. On the other hand, as we 
have shown in Fig. 5, it rapidly increases to b re i <~ 50 with the 
sensitivity to R <~ 26 in the uniform jet model. 

However we may have orphan afterglows due to an insuffi- 
cient gamma-ray sensitivity. Since the gamma-ray luminosity 
per unit solid angle rapidly increases with decreasing f? b s in the 
universal jet model, the detection of GRBs at large distances 
might be biased toward those with small O ^ S . If this is the case, 
we expect a much higher true GRB rate and large b K \ might be 
possible also in the universal jet picture. A test for this case 
is to see an anticorrelation between z and 6 \, s , and confirm or 
reject the 9 D \, S distribution obeying N(< # bs) oc (1 -cos0 o b s )/2, 
as predicted by the universal jet model, for a complete sample 
of GRBs within a redshift range. The distribution of 6 t> s can 
be observationally inferred from either luminosity function of 
GRBs with z measurements or jet break time obtained by after- 
glow light-curve fitting. However, the present sample of GRBs 
with known redshifts is too small and we must await future ob- 
servations. 

Another test possible by orphan afterglow observation is to 
examine the early behavior of orphan afterglow light curves. 
For # bs > #jet, in the uniform jet model we always expect no af- 
terglow flux at the earliest stage, and the afterglow flux should 
show gradual increase until the peak flux time, t p , given in §2. 
When b m \ » 1, the orphan afterglow events should be domi- 
nated by such cases, and hence we expect that the majority of 
orphan afterglows should show a slow rise at the beginning. 
On the other hand, the X-ray and optical emission should start 
immediately after the prompt burst when o b s < 9 1 in the uni- 
versal jet picture, since there is ejecta moving toward the ob- 
server. The initial slow rise of light curves is possible only when 
#obs > 7 , but such events should be relatively rare because of 
the large 9 1 inferred from observations and low absolute lumi- 

6 Rossi et al. (2002) have shown that the jet initial opening inferred in the un 
universal jet model. 



nosity for the large 9 \, s cases. The detection of such early rise 
of light curves may not be easy even in the case of the uniform 
jet model because the rising time scale is still smaller than the 
overall variability time scale T as can be seen in Figs. 1 and 2. 
If detected, however, it would argue for the uniform jet picture. 

6.3. Caveats of Our Predictions 

Here we describe several caveats of our prediction which 
should be kept in mind when one compares it to observed data. 

It is theoretically conceivable that the GRB central engine 
ejects not only the ultrarelativistic outflow that is responsible 
for GRBs, but also less relativistic matter with comparable to- 
tal energy. Prompt gamma-ray emission may be very dim or 
completely absent from such less relativistic components of 
ejecta, while the afterglow emission similar to those associ- 
ated with GRBs might be possible. Even if the ultrarelativis- 
tic component may be strongly collimated to produce beamed 
GRBs, more isotropic, less relativistic component could be as- 
sociated to most of GRBs. It is also possible that there are 
much greater number of events of dirty fireball or failed GRBs 
without prompt gamma-ray emission to any direction, than that 
of observed GRBs. The brightest class of core-collapse super- 
novae, called hypernovae (Iwamoto et al. 2000; Nakamura et 
al. 2001; Mazzali et al. 2002), might also be failed GRBs. 
In either case, the orphan afterglow rate can be increased sig- 
nificantly from our prediction. Even if an orphan afterglow is 
discovered, it is not an easy task to discriminate such another 
component of ejecta or failed GRBs from the pure effect of 
GRB jet collimation (Huang, Dai, & Lu 2002). 

One possible way of discrimination between an orphan af- 
terglow due solely to the viewing geometry from one due to 
a "dirty" fireball is provided by the afterglow decay rate. In 
the former case, the jet collimation should yield a decay slope 
a ~ -2, while a less collimated component or a failed GRB 
are expected to yield a ~ — 1. The slower rise of the light- 
curves of off-axis, collimated jets is also useful as discussed in 
the previous subsection (see also Huang et al. 2002). Stronger 
linear polarization is also expected for off-axis, collimated jets 
(Sari 1999; Ghisellini & Lazzati 1999; Granot et al. 2002). Fi- 
nally, the wavelength and sensitivity dependence of b m \ might 
be able to discriminate the above possibilities. The signature 
of off-axis afterglows would be the continuous increase of b K \ 
with the search sensitivity, though this trend might be mim- 
icked by failed GRBs if their event rate is continuously increas- 
ing with decreasing energy output to the dirty component. The 
constraint of b K \ < several by the past X-ray searches already 
indicates that the event rate of failed GRBs producing X-rays 
cannot be much higher than that of successful GRBs. On the 
other hand, the so-called X-ray rich GRBs might be a popula- 
tion between the failed and ordinary GRBs (e.g., Kippen et al. 
2002). 

We have assumed that the sky coverage of the gamma-ray 
emission is that corresponding to the jet opening angle obtained 
from afterglow modeling. However, this is not necessarily war- 
ranted. The causally connected angular scale during the prompt 
GRB phase is only 9 <~ T" 1 , which should be much smaller than 
the jet opening angle 0j et (Kumar & Piran 2000). Then it is pos- 
sible that the jet outflow, and thus its gamma-ray emission, is in- 
homogeneous on an angular scale much smaller than 0j et , while 
the afterglow emission arising later in the jet evolution is more 

jet model is in fact the observer's angular off-set relative to the jet axis in the 
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homogeneous. In this case the GRB rate is significantly under- 
estimated, as there could be a number of undetected GRBs due 
to some low-energy patches moving toward the observer, even 
if the observer is located within the jet opening angle. 

We have also assumed that the ten GRB afterglows used here 
are representative for all varieties of GRB afterglows. Although 
X-ray afterglows were observed for the majority of GRBs, op- 
tical afterglows were found for less than half of GRBs. Some 
afterglows were missed because of late follow-up observations, 
but there were also unusually dim optical afterglows (Fynbo et 
al. 2001; Lazzati, Covino & Ghisellini 2002). Several expla- 
nations can be considered including dust extinction and very 
large redshift. Recent studies of high redshift galaxies sug- 
gest that about half of stars are formed in very dusty galaxies, 
the fraction being much larger at high redshift than in the lo- 
cal universe (see, e.g., Totani & Takeuchi 2002; Ramirez-Ruiz, 
Trentham, & Blain 2002). The sample of the ten afterglows 
used here is clearly biased toward afterglows occurring in less 
dusty galaxies. Then our estimate of optical orphan afterglow 
rate may be overestimated by a factor of about two. When a 
GRB occurs in a molecular cloud or high density region with 
significant amount of dust, the dust along the direction of the 
jet might be destroyed by strong optical-UV flash (Waxman 
& Draine 2000) or early X-ray radiation (Fruchter, Krolik, & 
Rhoads 2001), making optical afterglows visible for on-axis 
observers. However, this effect should be small for off-axis 
observers, because the off-axis X-ray afterglow flux is much 
weaker when the light-curve peaks. This phenomenon may ef- 
fectively reduce the observed relative beaming factor b m \. Fi- 
nally, it should be noted that the prediction for radio orphan 
afterglows is not affected by dust extinction. 

7. CONCLUSIONS 

We presented a quantitative prediction for the detection rate 
of orphan GRB afterglows, based on one of the latest after- 
glow models that has been tested with a number of observed 
afterglows. We found that the orphan afterglow rate sensitively 
depends on afterglow model parameters, and a fairly large b K \ 
(>100) is possible for some types of GRBs, by an optical search 
with reasonable depth (R > 24). We derived our best-guess pre- 
diction of orphan afterglow rate, by taking a weighted mean of 
the ten sets of afterglow parameters that fit to ten well-observed 
afterglows. Although there are a number of effects or caveats 
which could significantly change our predictions, the predic- 



tion will be useful as "a baseline model" when we interpret the 
results of past and future surveys for extragalactic transient ob- 
jects. 

Our prediction is consistent with all the past surveys in X- 
rays, optical, and radio wavebands. Greiner et al. (2000) re- 
ported that there are about 10 possible candidates of orphan 
afterglows, which is interestingly very close to our expecta- 
tion. Although the SDSS transient reported by Vanden Berk 
et al. (2002) is very likely to be a radio-loud AGN, their search 
has already reached a meaningful sensitivity, since our expec- 
tation is about 0.2 afterglows for this search. A recent search 
by Levinson et al. (2002) found about 30 candidate radio after- 
glows by comparing the FIRST and NVSS surveys, while our 
model expects about two orphans. 

Detection of orphan afterglows seems not extremely diffi- 
cult in future surveys. Accumulation of data of advanced X- 
ray satellites such as Chandra and XMM-Newton might detect 
orphan afterglows whose rate is enhanced by a relative beam- 
ing factor of b K i <~ 40 compared with that for GRB-associated 
afterglows. OGLE III and ROTSE-III projects could detect a 
few and a few tens of orphan afterglows with b K \ ~ 3-4 in half 
a year at R <~ 20 and 19, respectively, providing nightly light 
curves. The southern SDSS observation (R ~ 23) could detect 
about 40 orphan afterglows with b m \ ~ 14 during five-year op- 
eration. Further deep optical surveys by Subaru/Suprime-Cam 
(R ~ 26) might detect an orphan afterglow with b K \ ~ 50. In 
the more distant future, GAIA could detect <~ 900 afterglows 
down to R ~ 20, and DMT could detect 500 and 2000 orphans 
at R = 24 and 29, respectively. ATA would detect about 200 
orphans with b m \ <~ 15 at ^0.3mJy (1.4GHz) in the radio band. 
Further effort for these searches is encouraged in the near fu- 
ture. Finally, we should make a point that a future GRB mission 
monitoring a significant part of all sky, like BATSE, is desired 
to check whether a candidate orphan afterglow is really an "or- 
phan". It is crucial especially for shallow searches of orphan 
afterglows with small relative beaming factor. 
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Table 1 



Expected Number of Afterglows in Various Surveys 



Survey Name 


Sensitivity 


Area [deg 2 ] 


A'exp 


N on 


N exp /N on 




OJogz 


1Q(log7-) 


Clog 7" 


(l/r)- 1 


X-ray Observation 


ROSAT a 


1 x 10~ 12 


9.4 x 10 4 


8.3 


2.8 


3.0 


0.93 


0.43 


1.5 


0.52 


0.81 


Optical Observations 


SDSS l b 


19 


1500 


0.20 


0.067 


3.0 


0.62 


0.50 


3.5 


0.68 


1.3 


SDSS 2 C 


19 


10000 


1.3 


0.45 


3.0 


0.62 


0.50 


3.5 


0.68 


1.3 


ROTSE-III 


19 


3.9 x 10 5 


53 


18 


3.0 


0.62 


0.50 


3.5 


0.68 


1.3 


OGLE III 


20 


5000 


2.2 


0.51 


4.3 


0.65 


0.48 


6.8 


0.62 


3.0 


GA1A 


20 


1.7 x 10 6 


720 


170 


4.3 


0.65 


0.48 


6.8 


0.62 


3.0 


Schaefer d 


21 


1500 


1.8 


0.29 


6.2 


0.69 


0.46 


12 


0.58 


5.7 


SDSS 3 e 


23 


4500 


36 


2.6 


14 


0.77 


0.42 


32 


0.51 


18 


Supernova' 


23 


200 


1.6 


0.11 


14 


0.77 


0.42 


32 


0.51 


18 


DMT I s 


24 


20000 


480 


21 


23 


0.81 


0.39 


58 


0.47 


36 


Subaru 


26 


5 


0.73 


0.013 


55 


0.93 


0.32 


150 


0.40 


97 


DMT 2 h 


29 


1000 


2100 


11 


190 


1.3 


0.23 


670 


0.32 


460 


Radio Observation 


FIRST/NVSS 1 


3.3 


5990 


2.0 


0.13 


16 


0.13 


0.32 


220 


0.39 


140 


ATA 


0.16 


5200 


210 


14 


15 


0.63 


0.32 


250 


0.38 


160 



Note. — Col. (2): Sensitivity in ^F^/(erg cirT 2 s~') at 1 keV for the X-ray band, in R magnitude for the optical band, and F„/mJy at 5GHz in the radio band. The 
sensitivity is corrected based on a typical afterglow spectrum when observation was made in a slightly different band. Col. (3): Area of a snapshot observation. When 
a survey is a consecutive monitoring longer than the typical afterglow time scale T, we converted the exposure (area X time) into an equivalent area of a snapshot 
observation using (l/T). (See text for detail.) Col. (4):The total number of all GRB afterglows, including orphans, expected to be detectable by a snapshot observation 
with the surface are shown in the third column. Col. (5): The same as Col. (4), but only for afterglows associated with prompt gamma-ray emission. Col. (6): The 
relative beaming factor b m \ = A^p/Non- Cols. (7, 8): The mean and la dispersion of logz. Cols. (9, 10): The mean and la dispersion of log7\ where T is the time 
duration over which afterglows are brighter than the sensitivities, in units of days. Col. (11): The mean of r _1 . 
"The ROSAT All Sky Survey (Greiner et al. 2000). 

b A search by Vanden Berk et al. (2002) using the first 1500deg 2 field of the SDSS data. 
c The same as SDSS 1, but after the completion of the SDSS project. 
d A search made by Schaefer (2002). 

" The number expected in the southern sky after the completion of the SDSS project. 
f The number expected in past supernova surveys. 
g The all-sky survey mode of the Dark Matter Telescope. 
'The deep probe mode of the Dark Matter Telescope. 

' The search made by Levinson et al. (2002) by using the FIRST and NVSS surveys. 
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Observation Time [day] 



FIG. 1 . — Light curves of off-axis GRB afterglows in X-ray band (IkeV) for the ten well-observed GRBs (name and the opening angle of the jet indicated in each 
panel). Note that the distance is assumed to be z = 1 for all GRBs, for comparison. The light curves are shown for different viewing angles from the center of the 
jet, as o b s = 1, 3, 5, 10, 20, and 30° for the solid, dotted, short-dashed, long-dashed, short-dot-dashed, and long-dot-dashed lines, respectively. 



12 




1 10 100 1 10 100 
Observation Time [day] 



FIG. 2. — The same as Fig. 1, but for the optical (R) band. 
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FIG. 3. — The same as Fig. 1, but for the radio (5GHz) band. 
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FIG. 4. — Expected number of orphan afterglows detectable by a snapshot observation per all sky, in the X-ray band (IkeV) as a function of sensitivity limit 
(upper panel). The solid line (N SX p) is for all afterglows including orphans, while the dashed line (Ntm) is for on-axis afterglows which are associated with observable 
prompt GRBs, i.e., those with 8 \, s < £?j et . In the lower panel, the ratio of N exp /N on , the mean values of redshift (10( logz ) , dashed line) and time duration (T in days) 
over which the flux is above a given sensitivity are shown. For the time duration, two different means of 10< logr > (short-dot-dashed) and {l/T)~ l (long-dot-dashed) 
are shown. 
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FIG. 5. — The same as Fig. 4, but for the optical (R) band. In addition, the expected number of supernovae of type la (dot-dashed line) and type II (dotted line) 
calculated by Woods & Loeb (1998) are also shown in the upper panel. 
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FIG. 7. — Contribution of each of the ten GRB afterglows to the detection rate N exp (shown by the thick solid line), in the three wavebands. The line markings for 
the ten afterglows are shown in the figure. 
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FIG. 8. — The optical light curves of orphan afterglows, for a fixed viewing angle 9 \, s = 20°, but various jet opening angles, 9j et . The solid, dotted, short- and 
long-dashed, and short- and long-dot-dashed lines are for 0j et = 1, 2, 3, 5, 10, and 20°, respectively. The redshift is assumed to be 1. Other parameters of the 
afterglow model are: £j et = 3 X 10 50 erg, n e x = 1cm -3 , e e = 0.05, tj = 3x 10~ 3 , and p = 2. 




FIG. 9. — The 0j et dependence of the relative beaming factor b Ie \, i.e., the ratio of rate for all afterglows, including orphans, to those associated with prompt GRBs. 
Different curves correspond to different search sensitivities in optical band, as indicated in the figure. Other parameters of the afterglow model are: £j et = 3 X 10 50 
erg, «ex = 1cm -3 , e e = 0.05, e B = 3 X 10" 3 , and p = 2. 



